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Qh, 1. Introduction 

j^. Heavy quarkonium systems have played a prominent role in our current understand- 

(-H \ ing of the Standard Model. Indeed, both charm and bottom quantum numbers were 

J^ ' discovered through heavy quarkonium systems, J/V^, the lightest vector charmo- 

• i-H , nium, and T(1S'7^, the lightest vector bottomonium, respectively. They have also 

^\ ' been important for our understanding of QCD, the sector of the Standard Model 

jrt ' which concerns strong interactions. Indeed, since the heavy quark masses m are 

larger than Aqcd, the typical hadronic scale, two important properties of QCD 
play a role in these systems, asymptotic freedom and confinement. On the one 
hand, asymptotic freedom explains the narrow width of the lower laying states^. 
On the other hand, they are the closest objects in nature to two ideal static color 
sources, whose energy behavior at large distances serves as an order parameter for 
confinement (in the absence of light quarks pi 

It was soon realized that, due to asymptotic freedom, for sufficiently heavy 
quark masses, heavy quarkonium systems should be similar to positronium and 
amenable for a weak coupling analysis^. Unfortunately, actual charm and bottom 
masses turned out not to be sufficiently heavy as to allow to explain the observed 
spectrum in the weak coupling regime^. However, they appeared to be heavy enough 
as to allow for a good phcnomcnological description of the spectrum by means of 
simple non-relativistic potential models (see Ref. Elfor a review). What to take as 
the potential was, and still is, the main input of such models. The question then 
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arose whether such potential could in principle be obtained from QCD if reliable 
non-perturbative techniques were at hand. Formulas were produced for it in terms 
of expectation values of Wilson loops (to be evaluated non-perturbatively) in a 
1/to expa nsion up to order 1/rn?, including spin dependent and velocity dependent 
terma^l^. However, when one loop results for the potential became available from 
direct QCD perturbative calculations^, it was realized that some of them were not 
correctly reproduced by the perturbative evaluation of the Wilson loop formulas. 

Inclusive heavy quarkonium decays to light particles were calculated using a 
factorization hypothesis. Namely, the short distance annihilation process was com- 
puted in QCD at quark level, and the long distance (non-perturbative) effects were 
taken into account by the wave function (or derivatives of it) at the origin, which 
was calculated using potential models or dropped from suitable ratios. However, 
again, when one loop results became available it was noticed that IR divergences 
appeared in some of the short distance calculations^. 

The understanding of the IR divergences was possible due to the introduction of 
Non-Relativistic QCD (NRQCD)ii4 It was shown that color octet operators, which 
are absent in potential models, were necessary to cancel the above IR divergences, 
and hence the factorization hypothesis used so far were wrong^^. The long distance 
part needs not only wave functions at the origin but also matrix elements of color 
octet operators, which were not computable in terms of potential models. 

Thus the lesson seemed to be that potential models cannot incorporate all rele- 
vant features of QCD for heavy quarkonium systems. However, a indication occurred 
that it may not be necessarily the case. If one recalculates the formulas for the QCD 
potential in terms of Wilson loops from NRQCD instead of directly from QCD, the 
discrepancies with the direct QCD calculation mentioned above disappeared, if the 
matching coefficients of NRQCD were calculated at one loop^. 

One of the aims of this brief review is to illustrate that suitable potential models 
can indeed be regarded as effective theories of QCD, and hence totally equivalent 
to it, in a very particular kinematical regime, and, as such, NRQCD color octet 
operators have a precise representation in them. This produces a number of model 
independent results for the inclusive decay widths to light particles and for the 
NRQCD matrix elements. The second aim is to illustrate that in the weak coupling 
regime, which corresponds to a different kinematical situation, potential models are 
not an effective theory of QCD. This regime is well understood and a number of 
higher order calculations are available. 

Before entering the issues above, let us mention that heavy quarkonium physics is 
experiencing a revival. Recently, new states have been discovered and new proce sses 
have been measured, some compatible with theoretical expectations^, others nol!^, 
which is triggering theoretical research. We refer the reader to R,ef.ll5lfor up-dates 
of the current status of the field, to Ref. Elfor an extensive theoretical review and 
to Ref. II 71 for a recent experimental account. 



Model Independent Results for Heavy Quarkonium 3 

2. Heavy quarkonium as a non-relativistic system 

A system is called non-relativistic if, in the center of mass frame, the typical three- 
momentum of a particle p is much smaller than its mass m. This implies that the 
non-relativistic energy E :~ y rn^ + p'^ —m ~ p^ /m is much smaller than p. Hence a 
hierarchy of scales exist m >> p >> E, which may be exploited in order to simplify 
calculations. In addition other scales may also be important depending on the par- 
ticular non-relativistic system. For heavy quarkonium, Aqcd is also important. In 
fact, it already enters in the definition of heavy quark, namely a quark whose mass 
fulfills TO >> Aqcd- Such a definition together with asymptotic freedom, which 
implies asim) << 1, suggests that heavy quarkonium, namely a heavy quark and 
a heavy antiquark (not necessarily of the same fiavor), is indeed a non-relativistic 
system. 

Rather than exploiting the inequalities m » p >> E and to >> Aqcd in every 
individual calculation, it is more convenient to built effective field theories (EFTs), 
which implement them at the Lagrangian level. This is the approach we will follow. 



3. Non-Relativistic QCD 

The Non-Relativistic QCD (NRQCD) Lagrangian has the following aspecll^ 



J „ D^ (T B [D,E] (T fDx,E] 

i^NRQCD = tp'<iDo+ \-CFg— \- CD g-z-^- + ics g — ^ 

Zm Zm 8m^ Hm'^ 
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" - I ^ ~ +JCS5 z — T 



+'''V'^"^2^,-''''^^ + '''''^ir+'''' Sml +-"^^ 



m'^ m'^ m'^ m'^ 

m^ m^ m^ to** 
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where 



08(1^0) = V^T-^XX^T-V , OseSi) = V'tT-o-xx^T-rrV' 

Oi(iFi)-Vt(_|H)x-x^(-§D)^ (2) 
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Oi(3Po) = iV't(-|D-^)xxt(-iD •^)V' 
Oi(3Pi) = i VU-f D X a)x ■ X^(-§ D X ,t)V' 

Oi(3p2) = Vt(_|H(VJ"))xxt(-i D (V^))V' . 

■0 is a Pauli spinor which annihilates a heavy quark and x a PauH spinor which 
creates a heavy antiquark. cf, cd, cs, /i, /s, etc. are matching coefScients which 
encode (non-analytic) contributions from (relativistic) energy scales of order m, and 
may have a factorization scale (/i) dependence. The NRQCD Lagrangian is obtained 
from QCD by integrating out energy fluctuations about the heavy quark mass and 
three-momenta higher than, or of the order of, m for the heavy quarks, and four 
momenta higher than, or of the order of m, in the gluo n fields . This can be done 
in perturbation theory in as(m/^ since as(m) << 1 (se a I I for an efficient way 
of doing such a calculation). Hence NRQCD is equivalent to QCD at any desired 
order in asijn) and 1/m. Note that the NRQCD Lagrangian is organized in inverse 
powers of to, which means that only the hierarchy m » Aqcd, P, E has been 
exploited. Hence, any dimensionful field in it does not have a definite size but may 
take the value of any of the remaining scales (Aqcd, P, E). In spite of this, a concrete 
velocity [v) counting was put forward in the original papers under the assumption 
that Aqcd ^ E —: rav"^ (then p ^ mv) which was useful to systematically organize 
calculations. As we will make clear in the following sections this is only one of 
the various counting possibilities that NRQCD admits, and may not be suitable 
for all heavy quarkonium states. In any counting, however, the scale dependence 
of the matching coefficients cancels against the scale dependence induced by UV 
divergences in NRQCD calculations, and, hence, each /x-dependence is eventually 
traded for one of the remaining dynamical scales (Aqcd, P, E). 

Note that the NRQCD Lagrangian is manifestly invariant under rotations but 
not under Lorentz transformations. The Lorentz symmetry is, however, non-linearly 
realized and provides constraints on some of the matching coefficients (for instance, 
cs — 2cp — l). These constraints were first uncovered in R.ef. 1201 using reparametriza- 
tion invariance. In Ref. |^it was shown that they follow form the Poincare algebra. 

The NRQCD Lagrangian contains non-hermitian terms due to imaginary pieces 
of the matching coefficients of the four quark operators (/i, /§,..) (see Ref. Hjfor a 
recent update). This is due to the fact that a heavy quark and a heavy antiquark of 
the same flavor may annihilate into hard giuons (of energy ~ m) , which have been 
integrated out. These non-hermitian pieces must be there in order to guarantee the 
equivalence of NRQCD and QCD. They contain crucial information about inclusive 
decay widths to light degrees of freedom. For instance, for P-wave states one obtains 
at leading order in the original velocity counting 
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+Im fsC''+'Ss){xQ{nJS)\OsCSo)\XQ{nJS))]. (3) 

where XQij^JS) stands for a heavy quarkonium P-wave state of principal quantum 
number n, total angular momentum J and spin S. This is to be compared with 
the potential model result, which is recovered by dropping the second term and 
identifying 



{XQ{nJS)\0,{^^+'Pj)\XQ{nJS)) = ^|i?(")'(0)p (4) 

zvr 

where R\-^{{y) is the wave function. The second term in (jSJ is however crucial in 
order to cancel the scale (/x) dependence of the matc hing: coefhcient of the first term 
at one loopl^. For instance, in the ^Pq case it reada^^ 



3D^ _ qr^„ / CU _/^„ \ _, ^0^^2 I 1 . ">5 



Im/i(^Po) = 3Cf \-f- - Cp\ ^a,{2mYll+ ^ 

/427 ^2 X 4 / 29 , 11 

+ (^-144)^-^+27"n"T-^°^'2;^ 



3 + -)C.+ 



(5) 



Let us also mention that the NRQCD formalism has also been used for the 
description of semi- inclusive decays (see Ref . |^ and references therein) and inclu- 
sive production (see Ref. EH and references therein) . We will not discuss these two 
applications here. 



4. Potential NRQCD 

NRQCD does not take advantage of the inequality p >> E.ln particular, it contains 
gluons of typical energy p, which cannot be produced in processes at the energy scale 
E. Simplifications should occur if one further integrates out degrees of freedom with 
energies larger than E, which leads to Potential NRQCD (pNRQCD)"^^. 

Unlike in NRQCD, the degrees of freedom, and hence the Lagrangian, of pN- 
RQCD depends on the interplay of Aqcd with p and E. We shall discuss two 
situations below: the weak coupling regime {k » E > Aqcd ) and the strong 
coupling regime {k > Aqcd >> E), where k is the typical momentum transfer, 
which, for low lying states, is of the order of p. 

^The /^-independent piece of this result slightly differs from the earlier calculations in Refs.l^ 
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4.1. Weak coupling regime 

li k >> E > Aqcd we can first integrate out energies ^ k . The EFT thus o btained 
is pNRQCD in the weak couphng regime. It has the following aspecll ^ l ^ ( I 



:pNRQCD = Tr {st {ido -h,)S + +0t {iDo - ho) o} + 

+ Tr <^ Qtr • gE S + H.c. + ^ + — ^ 



(6) 



where S = ^(R, r,t) and O — 0(R, r, i) are singlet and octet wave function fields 
respectively, R is the center of mass coordinate (whose dynamics is trivial at lower 
orders and has been neglected above), r is the relative coordinate, and hg and ho 
are quantum mechanical Haniiltonians 

hs = -—-Cf^ + ----^^-^(ifieSo)-2S^f,CSo)-hCS,))] + 
m r 2 m^ \ / 






TO V 2 r 



(7) 



2 m? 



4fsCSo)-2S^fsCSo)-f8CSi)) 



Tgj projects on states of spin S and total angular momentum J (see Ref. 28; for a 
precise definition). Only gluons and heavy quarks of energies smaller than k '^ 1/r 
are present in 0. However the only constraint on the three- momentum of the 
heavy quark is still that it must be smaller than m. The potentials in play 
the role of matching coefhcients. They can be calculated perturbatively in as{k) 
{as{k) << 1 since k >> Aqcd) and in the l/m expansion. Beyond tree level, 
this calculation produces UV divergences and also IR divergences if the smaller 
scales {E, Aqcd) are expanded. Once properly renormalized, the former cancel 
part of the scale dependence of the NRQCD matching coefhcients. The latter and 
the remaining scale dependences from NRQCD matching coefficients cancel with 
scale dependences induced by properly renormalized UV divergences in pNRQCD 
calculations (see Refs. l29ll3(Jl for illustrations in QED) . Note also that if we drop the 
octet field we recover a particular potential model. This is enough if we neglect non- 
perturbative contributions (meaning contributions for w hich th e scale Aqcd plays 
a role), and are only interested in corrections up to 0{af)^^^^. Beyond that order 
or if we wish to take into account non-perturbative contributions, the remaining 
gluon and octet fields are crucial. 
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Note also that, analogously to NRQCD, the pNRQCD Lagrangian is manifestly 
invariant under rotations, but not under Lorentz transformations. The constraints 
from the full Poincare algebra have been worked out in Ref. |2II 

Let us mention here that, in spite of the fact that the top quark decays due to 
the weak interactions before forming hadronic states, this regime is also relevant for 
the study of the top-antitop system near its production threshold""'. 

4.1.1. Spectrum 

The corrections to the spectr um up to order a^ had already been obtained before 
the introduction of pNRQClJ"^ ( and, in fact, making no use of NRQCD). The 
one loop potentials of the singlet field had been directly calculated from QCIJ"'. 
The two loop potential, which was also necessary at this order, was calculated 
using static heavy quark propagators^. NRQCD and pNRQCD just make the 
calculation simpler. Beyond that order or if one is interested in non-perturbative 
contributions due to the scale Aqcd all degrees of freedom of pNRQCD play a role 
and correct results cannot be obtained by just calculating potentials to a higher 
order''. In order to proceed further one has to specify the size of Aqcd with respect 
to E. If Aqcd ^ E, the leading non-perturbative effects are parameterized by 
non-local condensates and compete in size with the a^ perturbative corrections. 
If £^ >> Aqcd, one can carry out weak coupling calculations with the (ultrasoft) 
gluons in ©. The physical observables can then be organized in powers of as (at 
different scales) since p ~ /c '^ mas and E ~ mas^. The logarithmic contributions 
to the corrections at 0{a^) were calculated in Ref. "36" and the finite parts for the 
ground state in Ref. 1371 Not only that, the use of EFTs, in this case NRQCD and 
pNRQCD, allows to resum IR QCD logarithms. For heavy quarkonium systems this 
was first proposed in Ref. (12 within NRQCD, later addressed in a slightly different 
EFT framework called vNRQCE^, and implemented in the NRQCD-pNRQCD 
framework in Refs. I39I I4U I which produced the first correct NNLL resummations for 
the complete spectrunJ^ (see also^^ . NLL resummation for the hyperfine splitting 
have been obtained recentlys^. Non-perturbative contributions are parameterized 
in this case by local condensates^. 

4.1.2. Inclusive Decays 

The information on the parton subprocesses of inclusive decay widths to light par- 
ticles (light hadrons, photons or leptons) is encoded in the imaginary parts of the 
NRQCD matching coefficients. These are inherited in pNRQCD as imaginary parts 
of local potentials ((5(r) and derivatives of it), which eventually makes the decay 
width proportional to the wave function at the origin (or derivatives of it). Correc- 
tions up to Qfg to the wave function at the origin were known before the introduction 

''However, the leadings and next-to-leadin g^^ l non-perturbative contributions in the case E >> 
Aqcd were obtained before the introduction of pNRQCD. 
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of pNRQCE^^. In the case E » Aon n, the double logarithmic corrections at a: 



3 



were already obtained in this frameworlP^and, today, the single logarithmic correc- 
tions at that or der are also knowrl^. The resummation of logs at NLL has already 
been carried out . Semi-inclusive radiative decays have recently been addressed 

within this framework^. 



4.2. Strong coupling regime 

If fc > Aqcd >> E, the integration of all degrees of freedom with energies larger 
than E cannot be carried out in an expansion in ctg- In this case the degrees of free- 
dom of pNRQCD are the singlet field interacting with a potential and the pseudo- 
Goldstone bosons . If the latter are ignored, the form of the pNRQCD Lagrangian 
reduces to that of a potential model. 



-pNRQCD 



= Tr|s^(i9o-Ms| (8) 



where /i is a quantum mechanical Hamiltonian. Again, h is manifestly invariant 
under rotations, but not under Lorcntz transformations. The constraints of the full 
Poincare invariance have been discussed in Ref. 1471 

In the particular case k » Aqcd, we may integrate out first energies of the 
order of fc, which can be done in an expansion in as{k) and 1/m exactly in the same 
way as in the weak coupling regime. We are thus lead to the same Lagrangian ©, 
which may be renamed as pNRQCD' since it is not our final EFT yet. We still have 
to integrate out energies ~ Aqcd- This cannot be done in perturbation theory in 
tts anymore but one can use the fact that k, p » Aqcd >> E. If one assumes 
that the octet field develops a gap ^ Aqcd , it can be integrated out and we are left 
with the singlet field only. Namely we recover the degrees of freedom of a potential 
model as explicitely shown in ^. 

In the general case fc '^ Aqcd , the integration of energies of order fc cannot be 
done perturbatively in ag. If one assumes that the potentials are analytic in 1/m, 

V^ V, V-, Va 

m m ra'^ m, 

one can obtain them by matching IJSJ to NRQCD in the l/m expansion. Then one 
obtains the non-pert urb atiye potentials in terms of expectation values of operator 
insertions in Wilson loopa^. In this way one is able to reproduce and correct earlier 
results^. In particular, it was noticed in this approach that the 1/m potential had 
been missed before. 



/ .g2 

V] (r) — lim — - — 
^ ' T-,oo \ AT 



T/2 j-T/2 

dt dt'\t-t'\ 
-T/2 J-T/2 




-1 


((E(i) . E(f')))n - 


- ((E(t)»n . ((E(t')))n 




. 






(10) 
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where D stands for a rectangle of size T x r. ((• ■ • ))□ means the expectation value 
of the depicted fields joined by Wilson lines along the rectangle, divided by the 
expectation value of the Wilson loop in the same rectangle. 

Let us make a parenthesis here and exemplify how the mismatch between the 
earlier Wilson loop approach and the explicit QCD one loop calculations mentioned 
in the introduction is resolved in the present formalism. Consider, for instance, one 
of the terms in V2 contributing to the hyperfine splitting. 



V^l''\r) = ^i hm / dt ((gBi(i) • .gB2(0)))n - 4(4. + d.vCf) S^^^^r) 



(11) 



3 T- 



where the dsv and d^ are suitable combinations of the fii^^^Ss) and fsC^^^^Ss) 
matching coefficients of the four fermion operators in (1)^. In the earlier Wilson 
loop approach one would obtain the same expression with cf = 1 and dsv = dvv = 0, 
namely the short distance contributions coming from loops and virtual annihilation 
processes from scales of the order of m were missing. If one calculates the chromo- 
magnetic correlator at one loop one finds a contribution proportional to a^ log(r/x) 
which adds to a a^ log(?7i//i) contribution in d^v + d^C f producing the full QCD 
ttg log(r?Ti) contribution. 

Let us remark that the short distance behavior of these potentials can be calcu- 
lated in perturbation theory in as{r) and hence they must coincide with the ones 
in hs of the weak coupling regime (7). Therefore, they become increasingly singular 
at short distances as we go further in the 1/m expansion. Hence the Hamiltonian 
h is not well defined in standard quantum mechanics. In order to make sense of it 
we must understand this Hamiltonian as an EFT. As such we should regulate it, 
establish a counting, and treat the subleading pieces as perturbations. The scale 
dependence induced by the regularization should cancel exactly with that in the 
NRQCD matching coefficients, much in the same way as it was observed in Rcf. 1301 
for QED. 

The power counting in h depends on the typical value of p and r in the concrete 
bound state we wish to analyze, and hence a simple power counting cannot be fixed 
a priori. Only a few statements can be made in general. The kinetic term and Vb 
must always be assigned the same size {mv^ since p ~ mv) and taken as leading 
order. Although Vi is suppressed by a^ and hence order mv in the weak coupling 
regime, it may in our case also be leading order, since in the strong coupling regime 
as ~ 1 and dimensional counting allows for a size Vi ^ (mv)'^. The terms in V2 
are at most order mv'^ (although in the weak coupling regime are order mv'^ due to 
extra ctg suppressions) and hence they can be treated as perturbations. 

One should be aware that on general grounds the form of the potentials in (9) 
is not unique. Unitary transformations are allowed in quantum mechanics which 
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change the aspect of the Hamihonian but do not change the physics. In EFT one 
should better stick to transformations which preserve the counting. Even withm 
those one can reshuffle contributions from one term to another in the potentials^. 
In general different ways of performing the matching from NRQCD (or directly from 
QCD) lead to different forms of the potential related by unitary transformations. 
For instance, in the weak coupling regime, matching on-shell matrix elements rather 
than using the 1/m expansion, or matching in the Feynman gauge rather than in 
the Coulomb gauge produces different forms of the potential. 

Recently, it has been shown that contributions to the potentials which are non- 
analytic in 1/m exists. A procedure to compute the ones due to the three-momentum 
scale y^toAqcd was put forward in Ref. 0^1 They give rise to subleading contribu- 
tions with respect to the 1/m potentials. All these potentials (analytic and non- 
analytic) can be evaluated on the lattice . Further non-analytic terms may appear 
due to the three- momentum scale mas when this scale is much larger than Aqcd, 
which have not been taken into account so far. 



4.2.1. Spectrum 

Once the non-perturbative potentials ar e obtai ned from a lattice calculation (or by 
means of other non-perturbative method j^ ' ' I ' ' ' *), one may think that the Schrodinger 
equation can be solved and the spectrum obtainecp^in total analogy with potential 
models". A fully consistent calculation, however, requires the lattice calculation of 
the potentials to be translated to MS scheme or similar, in order to match the avail- 
able NRQCD matching coefficients, or vice- versa. Furthermore, because of the same 
reason, one has to use (or to translate to) MS scheme the calculations in quantum 
mechanics perturbation theory. The advantage is that now one has a counting and a 
well-defined procedure which allows, at least in principle, to systematically improve 
the calculation by adding higher order potentials and by going to higher order in 
quantum mechanics perturbation theory. The contribution to the spectrum from 
potentials which are non-analytic in 1/m turns out to be very suppressed. The role 
of pseudo-Goldstone modes has not been addressed in this framework so far. 

4.2.2. Inclusive decays 

As in the weak coupling regime, the imaginary parts of the NRQCD matrix elements 
are inherited in local terms of the pNRQCD Lagrangian. In particular, the local 
terms in hg (7) also exists in h. More important the color octet operators of NRQCD 
also have a representation in terms of local potentials. For instance, if we restrict 
ourselves to the P-wave contributions we have 

luiVi = CAT's'jV'5^^\v)Vnm fiC^+^ Pj) (12) 

-H^f3V,S(3)(r)V UlTnhCSo) - 2 S^ {ixn^CSo) - Im/sC^^i)) J + • • ■ 
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Then the decay width of P-wave states to hght hadrons at leading order now reads'^ 



TixQinJS) ^ LH) 



3Im hC'+'Pj) + l^lm fs{''+'Ss)£3 

oCa 



(13) 



where £3 is a non perturbative parameter defined as 



1 1"^ 
£3 = —/ dtt^(gE{t)-gE{0)), (14) 

^^c Jo 

By comparing with (Q , one may rephrase (|13|l as 

{XQ{nJS)\OsCS,)\xQ{nJS)) = -^ '^"^ V ^3- (15) 

Namely, one is able to obtain NRQCD color-octet matrix elements in terms of 
(derivatives of) wave functions at the origin, which contain all flavor and princi- 
pal quantum number dependence, plus extra universal (depending on Aqcd only) 
non-perturbative parameters. One can check in perturbation theory that the scale 
dependence of £3 cancels exactly the scale dependence of Im fi{^^'^^Pj) (5). The 
unknown non-perturbative parameters together with the wave functions at the ori- 
gin may drop from suitable ratios. They may also be extracted from data. This 
allowed to put forward a pr ediction for bottomonium states in terms of data ex- 
tracted from charmoniurrl^, which turned o ut t o be in reasonable agreement with 
the experimental results when they came oulP^. 

If all the states below threshold for bottomonium and charmonium were in the 
strong coupling regime, and if we restrict ourself to potentials which are analytic 
in 1/m, one would o btai n a reduction of unknown NRQCD matrix elements by 
roughly a factor of two^^lat order 1/m^ (i.e. LO for P-wave states and NLO for S- 
wave states). Non-analytic terms ^ ->/toAqcd give rise to subleading contributions 
(provided that ^ttiAqcd >> rnus ( -y/m- Aqcd ) ) , which, however, may be of the 
same order as analytic corrections to the leading analytic contributiona^. They 
would slightly increase the number of unknown matrix elements. 

5. Discussion 

So far we have discussed a theoretical framework with almost no reference to any 
heavy quarkonium state in nature. If we wish to apply it to a concrete state we 
have to first figure out whether this state belongs to the weak or strong coupling 
regime, if any. This is not easy to establish a priori, since what scale plays the 
role of Aqcd or the typical value of k or even E, cannot be extracted directly 
from the experimental observables. One may try the weak coupling regime first and 
check: (i) if the expansion in ctg shows good convergence and (ii) if the leading non- 
perturbative effects are small. These appear to be fulfilled by the T(15') system. 
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and to a lesser extend for the B^ and J/tp '^. The constrain (ii) is very restrictive 
since non-perturbative effects in the weak couphng regime grow with a large power 
of the principal quantum number'^. Hence, most likely any excited state does not 
belong to the weak coupling regime. This does not mean that they belong to the 
strong coupling regime, since there is also the possibility that pNRQCD does not 
apply to them . Indeed, the constrain k ^ Aqcd >> E, does not allow states 
with k '^ E ^ Aqcd ■ If we accept Heavy Quark Effective Theory counting rulea"^, 
these are states close to, or beyond, the heavy-light pair production threshold. 
Hence neither the weak nor the strong coupling regime are in principle applicable 
to these states. In order to make things concrete, let us advance what we believe to 
be a reasonable (although not entirely conservative) assignment. For the T system, 
T(IS') and rjb^lS) belong to the weak coupling regime, whereas the remaining states 
below heavy-light pair production threshold may well be considered in the strong 
coupling regime. For the ip system, J/ip and r]c{lS) seem to be in the border line 
between the weak and strong coupling regimen, the lowest lying P-wave states in the 
strong coupling regime, and the remaining states (including tp{2S) and ric{2S)) are 
either too close or beyond the heavy-light pair production threshold, so most likely 
pNRQCD is not applicable. For the Be system, the pseudoscalar and vector ground 
states may well be in the weak coupling regime^ ' , the lowest lying P-wave states 
and the first radial excitation of the S'-wave states in the strong coupling regime, 
whereas for the remaining states pNRQCD would not be applicable. 

6. Conclusions 

EFTs techniques allow to exploit efficiently the various hierarchies of scales appear- 
ing in heavy quarkonium. They simplify and systematize earlier approaches ( for 
instance, potential models) and are powerful enough to put forward new results, 
which stem from QCD under well controlled approximations. This is so both in 
the weak and the strong coupling regime. In the weak coupling regime a number 
of explicit calculations have been carried out at higher orders in as . In the strong 
coupling regime non-trivial results for the NRQCD matrix elements have been ob- 
tained. The phenomenological consequences of these results have not been fully 
exploited yet. 
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